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We revisited a model for charmonium hybrid meson with a magnetic gluon [Yu. S. Kalashnikova 
and A. V. Nefediev, Phys. Rev. D 77, 054025 (2008)] and improved the numerical calculations. These 
improvements support the hybrid meson interpretation of X(4260). Within the same model, we 
computed the hybrid meson mass with an electric gluon which is resolved to be lighter. Relativistic 
effects and coupling channels decreased also the mass. 
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I. INTRODUCTION 

Quantum Chromodynamics (QCD) is widely believed 
to be the theory of strong interactions. In this non 
Abelian theory, the gauge bosons, the gluons, carry a 
color charge. We are therefore led to observable color sin- 
glet configurations made of quarks but also gluons. Glue- 
balls, bound states of only gluons, and hybrid mesons 
are consequences and predictions of QCD. Besides glue- 
balls, hybrid mesons deserve much interest. The gluonic 
excitation leads to low-lying states with quantum num- 
bers not allowed for usual mesons (Note that the low- 
lying gluon states are non exotic but glueballs with ex- 
otic quantum number also exits and are called oddballs, 
for a review see Ref. [![). Their observation would be 
another great confirmation of QCD. However, the detec- 
tion of exotic hadrons remains a challenging task (for a 
review see Ref. |2(). 

The properties of hybrid mesons were investigated in 
various approaches, see for examples Ref. [H, 0] , and in 
lattice QCD Q. The former have two common inter- 
pretations. In the first one, the flux tube linking the 
quark to the antiquark is in an excited state, allowing 
quantum numbers that cannot be reached by the usual 
quark-antiquark picture. In the other scenario, the ex- 
citation is modelled by a constituent gluon, leading to 
a three-body system. We know that both picture are 
closely linked: The constituent gluon creates an equiva- 
lent potential compatible with the energy of an excited 
string [|,0|. 

All these approaches rely on quasiparticle interpreta- 
tion. For heavy particles, the spectrum and decay prop- 
erties are extracted from a Hamiltonian formalism. One 
has then to resort to a numerical procedure to find eigen- 
values and eigenfunctions. For two-body systems, like 
heavy charmonia, different methods exist. The most ef- 
ficient one is the so called Lagrange mesh method. This 
technique allows one to compute energies and wave func- 
tions straightforwardly since the method is not varia- 
tional 1 81. Matrix elements for semi-relativistic kinetic 
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energy are easily computed. However, it does not admit 
a three-body generalization in the semi-relativistic case. 
For three-body systems, like hybrid mesons, one gener- 
ally computes the mass within a finite dimensional basis. 
We expand the unknown wave function on a set of trial 
functions forming a basis of the Hilbert space in the limit 
where the dimension of the basis goes to infinity. 

In practice, we only deal with a finite number of basis 
functions. Fortunately, we know the approximation of 
the mass to be always an upper of the true eigenvalue, 
allowing to minimize the mass with respect to parameters 
19!]. The question is then: How many basis functions 
should we consider to have an acceptable accuracy ? In 
this work, we answer to this question for the three-body 
description of hybrid mesons. 

Recently, Kalashnikova and Nefediev investigated the 
spectrum of the charmonium hybrid ccg within the con- 
stituent gluon model [lfj. They derived a Hamiltonian 
thanks to the field correlator method and introduced ein- 
bein (or auxiliary) fields to deal with relativistic kinemat- 
ics. The resulting Hamiltonian was diagonalized with one 
trial function that was taken as a Gaussian depending of 
a hyperradius. Finally, the authors computed correction 
coming from the string, self-energy and spin-dependent 
operators. The authors studied the case of the lowest 
ccg hybrid meson in with the gluon quantum numbers 
are £ g = 1 and j = 1 (magnetic gluon) . In this work, we 
investigate how accurate are those approximations for 
the wave function (truncation of the basis) and for the 
Hamiltonian (introduction of auxiliary fields) and greatly 
improve the accuracy on the mass thanks to correlated 
Gaussian functions. 

In Sec. [Ill we present the model of Ref. 1QJ and recall 
its main properties. The authors used a particular ansatz 
for the wave wave function which is described in Sec. IIIII 
We discuss some improvements for the wave function in 
Sec. IIV1 In this section, which consist in the main part 
of this work, we perform numerical calculations with var- 
ious approximations for the model and make some com- 
parisons. Having identified a good approximation for the 
wave function, we test in Sec. [V] the stability of the error 
under parameter evolution. Finally, we draw our conclu- 
sion in Sec. IVI] 
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II. HYBRID MESONS AS qqg SYSTEMS 

The starting point for the hybrid model of is a 
Hamiltonian derived from QCD thanks to the field cor- 
relator method. For a hybrid meson seen as a qqg system, 
this Hamiltonian reads 



H = H 



Vc, 



where 



H 



P q + 
<r\x„ 



Pa + m 2 g + 



1 

a\xr. 



V c 



3a s 



3a s 



2 1 "Eg "^q\ 



+ 



6 1 Xq Xq | 



(i) 



(2) 



(3) 



The first terms of Hq represent the kinetic energy of 
two quarks of the same mass and of the massless con- 
stituent gluon. The long-range confining force consists 
in two strings, two fundamental flux tubes of energy 
density c, linking the gluon to the quarks. Such a con- 
fining potential is in a gree ment with other phenomeno- 
logical approaches flll.Tl2l as well as with recent lattice 
QCD computations [l3|. For other three-body systems 
like baryons, the confinement is generally assumed to 
be a Y-shape. In the case of the hybrid, however, un- 
der Casimir scaling hypothesis for the string tension, the 
V-configuration is resolved to be energetically favorable 
The Torricelli point of the Y-shape merges with 
the gluon position. Imposing a Y-shape would obviously 
raise the whole spectrum. The Casimir scaling hypothe- 
sis is supported by recent lattice studies, see for instance 
ref. [13 . 

The confining potential is augmented by the short- 
range Coulomb potential Vc, where a s is the strong 
coupling constant, arising from one gluon exchange ef- 
fects between the constituent particles. In particular, 
the Coulomb potential between the quark and the anti- 
quark is repulsive since the qq pair is in a color octet. 
We refer the reader to Ref. [l5| for a detailed discus- 
sion about the short-range interaction potentials in QCD. 
Spin-dependent short-range interactions could be added 
in perturbation as it is done in Ref. [loj . However, it is 
enough for our purpose to consider only the dominant 
order Hamiltonian |T]). 

To be complete, we recall that the parity and charge 
conjugation of a qqg system are given respectively by 

p = c = (-l) 1+f 94 +s <«, (4) 

where i q g and i g are the orbital angular momenta of the 
qq pair and of the gluon, and where S qq - is the intrin- 
sic spin of the qq pair. Compared to usual mesons, ex- 
tra phases give rise to exotic quantum numbers such as 
1 h for example. Two cases have to be distinguished 
following the total spin j of the gluon: The magnetic 
gluon, for which £ g — j, and the electric gluon for which 



£ g = j ± 1 Although it may gain a dynamical 

mass induced by confinig forces, the gluon is a massless 
particle that remains transverse, with only two polari- 
sations P, [l6j]. Therefore, the minimal value for j is 1. 
Exotic quantum numbers require at least one P-wave. 
In the following, we shall refer to magnetic (electric), the 
qqg states with l g — 1 (£ g — 0). 

Solving the eigenequation associated to a three-body 
semirelativistic Hamiltonian such as {l} is a difficult nu- 
merical problem. The intrinsic complexity of three-body 
systems obviously comes into play, but another problem 
is the determination of the matrix elements for kinetic op- 
erators of the form \J p 2 + to 2 , which are not commonly 
found in quantum mechanics. This last difficulty can be 
avoided by introducing auxiliary (or einbein) fields to get 
rid of the square roots appearing in Hamiltonian (TTJ) . One 
obtains 



H(fJ. q ,flg,flg) 



fig + Hg + fig Pg+m^ q p Z g 



2fl q 



2flg 



a\x a - x„ + a\x a - x 



2/J.g 



(5) 



where fi q , fi q , and fi g are the so-called auxiliary fields. 
Being formally defined as operators, they can be elimi- 
nated through the following equations 



S^H(H q ,H q ,Hg) 





= = 






= = 


fa 


Mg=£g 


= = 


> fa 



Pa 



m„ 



P% 



(G) 
(7) 
(8) 



It is then readily checked that H(p, q , fi qi fl g ) = Hq; both 
Hamiltonians are equivalent up to the elimination of 
the auxiliary fields as operators. However, the calcu- 
lations are considerably simplified if one considers them 
are c- numbers variational parameters. The eigenvalues of 
the spinless Hamiltonian ([5|) are more easily found since 
only nonrelativistic kinetic operators are present. These 
eigenvalues, denoted as Eo(n q , \i q , \i g \ are finally mini- 
mized with respect to the einbeins. The optimal values 
of the auxiliary fields seen as variational parameters are 
logically close to the average values of the correspond- 
ing operators [13, [Hj]. For example, the optimal value 
of fig, denoted as /i s o, is such that /i s o « (As)- It can 
be interpreted as a dynamical gluon mass. The same 
arguments hold for the other auxiliary fields. An impor- 
tant point, that has been shown in Ref. [19j |. is that the 
eigenvalues of Hamiltonian ([5]) are upper bounds of the 
eigenvalues of Hamiltonian (fT]) . The more auxiliary fields 
are introduced, the less this bound is accurate. In partic- 
ular, the accuracy of the auxiliary field method decreases 
when light or massless particles are present in the system 
under study. 

In this work we focus on the case where the quark and 
the antiquark have the same mass, i.e. m q = m q = to. 
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Then, by symmetry, fi q — /iq = p and, using the Jacobi 
coordinates 



A Xg Xq-i P — Xg 



(9) 



Hamiltonian |(SJ| becomes 

U {P, Ms ) =P + it H 1 1" 

2 p p 2q> 





A 




A 


a 


2~ P 


+ cr 


2 + " 



(10) 



where pa and p p are the momenta associated to A and p 
respectively, and where 



_ 2flflg 

2H + fig ' 

The center of mass, defined as 

pX q pXq -\- PgXg 



R 



2/1 + fig 



(11) 



(12) 



is decoupled and its conjugate momentum is set equal to 
zero since we work in the rest frame of the system. 



III. HYPERSPHERICAL FORMALISM 

Eigenvalues of Hamiltonian H(p,p g ) with auxiliary 
fields (|10p in the case of peg systems have been numer- 
ically computed in Ref. [1 01 ] by using the hyperspheri- 
cal formalism. We recall in this section the procedure 
that has been used in this last reference. The authors of 
[l(| state that since the qq pair is heavy, the assumption 
Iqq = can be made. Moreover, the particular case of 
a magnetic gluon with i g — 1 is considered in this refer- 
ence. For this particular choice, the spin of the quark pair 
Sqq is a good quantum number and leads to the following 
states: 1~ for S q9 = and (0,1, 2)"+ for S q g = 1. All 
these states are degenerate since H(p, p g ) does not have 
any spin-dependent term. This leads to wave functions 
of the form 



with the hyperradius 



0,£ g = l) 1 ® \Sqq) 



X 2 = — A 2 



2^lg o 



2/i + /ig 



(13) 



(14) 



\iqq,i g ) is a shorthand notation for the coupling 
[yi qjj (X) Cg) ye g (p)Y, with the solid spherical harmonics 
ye g m(p) — P la Y£ g m(p)- The following states can conse- 
quently be described within this approach 



|1 — ) = \x(X))<2> [|0, l) 1 c» |0) 

= |x(x)>® ki) 1 ®!!) 



(15) 
(16) 



with J = {0, 1, 2}. The trial wave function is chosen to 
be a Gaussian depending on one variational parameter 
(3, i.e. 



xiX))®^!) 1 =exp(--f3 2 X 2 )y lm (p). 



1 



(17) 

This solid harmonic, ^im(p) = pYi m {p), determines the 
angular momentum in order to treat exotic mesons with a 
P-wave gluon and a S'-wave qq pair. The spin function for 
the quark pair is irrelevant since our Hamiltonian is spin- 
independent. We will therefore get the same mass for the 
the four states ^ 1 , 1 '", 2 ^ under consideration. 

A convenient method to fit the parameters is to repro- 
duce the charmonium spectrum within the same flux tube 
model, that is with the Hamiltonian H cS = 2^/p 2 + m 2 + 
ar — 4a s /3r. The parameters obtained in this way are 
compatible with typical values used in potential mod- 
els 0: 



a = 0.16 GeV 2 , a, = 0.55, 



= 1.48 GeV. (18) 



We, as the authors of ref. [Hj, then implicitly assumed 
the same value for the parameters for both systems (char- 
monium and hybrid charmonium). It is also natural to 
keep the same parameters as in ref. [Toj since the present 
paper focuses on the improvement of the numerical reso- 
lution. In Sec. El we vary the parameters and investigate 
this influence of the one Gaussian approximation. 

A numerical resolution of H (/x, [i g ) with the ansatz (|17j) 
leads to [l(| 

Mq = 4.573 GeV, p, = 1.598 GeV, n g0 = 1.085 GeV, 

(19) 

for the ground state. This mass corresponds in a first ap- 
proximation to cc hybrid mesons with quantum numbers 
1— and (0, 1,2)-+. 

Spin-dependent corrections have also been computed 
to be equal to floj 



AM(0 _+ ) = -321 MeV, 

AM(r+) = -253 MeV, 

AM(r _ ) = -176 MeV, 

AM(2"+) = -116 MeV. 



(20a) 
(20b) 
(20c) 
(20d) 



The splittings induced by the spin-dependent operators 
in perturbation theory are clearly not negligible. How- 
ever, we know from the Rayleigh-Ritz method @ that 
the mass Mq = 4.573 GeV can only be an upper bound 
of the true eigenvalue of H(fi,p g ). Indeed, when trun- 
cating the basis with a finite number of basis function, 
one restrict the operator to the subspace spanned by the 
trial functions. It is then natural to investigate how far is 
Mq from the true eigenvalue and compare the difference 
with respect to the additional corrections (f2"0"l) . 
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IV. CORRELATED GAUSSIAN BASIS 

The real eigenvalue M of a Schrodinger equation 
H\$>) = Af |$) is generally found by expanding the wave 
function $(A, p) = (p, A|$) in a basis 

N 

<&(A,p) = 5>^(A,p). (21) 

i 

The real eigenvalue is reached in the limit N — > oo. Since 
for fixed AT, is an upper bound, we can then min- 

imize the mass with respect to the parameters on which 
depend cp % . A appropriate ansatz for tp l will lead to an 
accurate mass for low values of N. 

For the trial functions, we use a generalisation of the 
hyperspherical function {XT]) 

(/(A, p) = exp (-a 4 A 2 - bip 2 - 2c;A • p) yLAi(d l X+e l p). 

(22) 

The wave function used in [lfj is recovered for the par- 
ticular values N = 1, c\ = di = and the condition ([M]) 
on the variational parameters ai et b±. The correlated 
Gaussian (j2"2"|) shared many advantages 9]. They admit 
a generalisation for an arbitrary number of interacting 
particles thanks to matrix notation; matrix elements for 
usual power laws potential are expressed in close forms; 
the Fourier transform of a Gaussian is also a Gaussian. 
In particular, we have shown that matrix elements for 
semirelativistic kinematics \/ p 2 + m 2 can be easily de- 
rived m eh. 

This form is very convenient for hybrid mesons since 
we deal easily with angular excitations for the gluon but 
also for the quarks. Imposing some restriction on the pa- 
rameters we get two ansatz for the lowest P-wave hybrid 
mesons 

<p A (\, p) = exp (-a\ 2 - bp 2 ) y lm (p), (23) 
ip B (X,p) = exp (-a\ 2 - bp 2 ) y lm {\). (24) 

The function lpa(\, p) has the angular part \£ q g — 0, l g = 
l) 1 and correspond to a hybrid meson with a magnetic 
gluon, and the function y>s(A, p) has the angular behav- 
ior \£ q g — \,£ g — Q) 1 and corresponds to a hybrid meson 
with an electric gluon. 



and b H — 0.156 GeV. The relation (fl4|) with the param- 
eters (1181) implies a nearly symmetric Gaussian shape, 
i.e 

% = = 1.0138. (25) 

Here we perform a numerical optimisation on the two 
parameters of our single wave function </?4_[A,p). We 
expect to find a mass slightly lower than in . Indeed, 
the resulting mass for the four states h , l -- , 1 h , 2 h 
in this approximation reads 

= 4.462 GeV. (26) 

We gain 109 MeV by simply choosing a more general 
ansatz (with two variational parameters instead of one 
parameter) for the trial wave function. The shift is of 
order of the spin-splittings (|20|) . Let us look at the wave 
in order to understand this difference. The optimisation 
procedure leads to oi = 0.083 GeV and b x = 0.204 GeV. 
Their ratio strongly differs from (j25|) . As an illustra- 
tion, we plot the spatial part of the wave functions in 
Figs Q] and [2] (remind the factor p in the solid spherical 
harmonic) . 

Wave function given by the hypersherical formalism 




FIG. 1: Wave function (|17p for a magnetic gluon hybrid me- 
son. 



A. P-wave hybrids with magnetic gluon 

Let us now see how changes the mass when relaxing 
the condition (|14j) on a\ and b\. In order to have a rel- 
evant comparison point, we consider PT(/i,/^ s ) with the 
parameters given in (fT5|) and (fil?)) . If we use the hy- 
perspherical formalism with only one trial function (and 
hence only one variational parameter) in our basis, we 
find a mass of Mn = 4.573 GeV (obviously the same 
value as in Ref. [l(|)- The optimisation procedure leads 
to/3 = 0.62 GeV 1 / 2 which correspond to a H = 0.154 GeV 



We can go even further by adding more functions in the 
basis. With a sufficient number of functions we should 
converge to the real eigenvalue of our operator H (/z, p, g ). 
Convergence is obtained for a small number of Gaussians. 
With TV = 10 functions, the numerical procedure is ac- 
curate up to 1 MeV which is clearly enough. The wave 
functions is then a sum of 10 Gaussian trial functions. 
The coefficients and parameters of each functions can be 
read in Table HI The resulting mass 

M^ 10) = 4.445 GeV (27) 
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Wave function with N = 1 function in the basis 



0.04. 




FIG. 2: Wave function (I23[l for a magnetic gluon hybrid me- 
son. 



TABLE I: Wave function coefficients for magnetic gluon ([23}. 



i 


a, 10" 3 




bi 


1 


26.183 


0.082811 


0.203578 


2 


116.410 


0.207870 


0.697459 


3 


-6.902 


0.435376 


0.232531 


4 


-3.242 


2.429518 


0.720520 


5 


-4.284 


0.050134 


0.088911 


6 


-11.033 


0.079373 


1.083885 


7 


46.785 


0.153637 


0.383596 


8 


8.5083 


0.053970 


0.097767 


9 


-0.505 


0.165116 


0.078640 


10 


-98.834 


0.208914 


0.571826 



BaBar collaboration (22J. But they mentioned also that 
their mass M + AM(1 — ) = 4.397 GeV was higher than 
the candidate's mass. It is worth mentioning that the ex- 
act value of the Hamiltonian H(fi,^, g ) used in is 130 
MeV below the approximative value computed with the 
trial function (JTTJ) . Assuming that spin-dependent cor- 
rections do not change much with the change of the wave 
function, we get a mass compatible with the X(4260). 

Wave function with N = 10 functions in the basis 




FIG. 3: Wave function ((22j) with 10 Gaussian functions J23J) 
for a magnetic gluon hybrid meson. 



B. P-wave hybrids with electric gluon 



is somewhat lower. The gain when increasing the ba- 
sis is not as strong as the first gain obtained with one 
trial Gaussian ([23)) with two variational parameters. We 
can understand this effect as follows: The Gaussian wave 
function ([2"3")) has a good overlap with the exact wave 
function of the three-body system. However, when we 
restrict the trial function to be a Gaussian of the only 
hyperradius (|17p . the overlap with the wave function is 
worst since the hyperradius imposes a relation relation 
between the two coefficients of Jacobi variables which is 
not realised in the optimal case. The wave functions with 
N = 10 functions is shown in Fig. [3J The relative error 
between the first trial function (JTTJ) and the eigenfunc- 
tion is 33% and the relative error (i.e. the ratio between 
the integral of the difference squared and the integral of 
the function with N = 10 squared) with the single Gaus- 
sian approximation is 18%. Nevertheless, the mass found 
with one hyperspheric Gaussian is less than 3% above the 
real value. 

In their paper, the authors of [l 0j | concluded favorably 
in a hybrid meson interpretation of the vector meson can- 
didate V(4260). This new state was observed by the 



In the previous section, we focused on the magnetic 
gluon hybrids. We followed the approach of [l(| where 
the authors assumed a S'-wave quark-antiquark pair. 
Since, the quark-antiquark pair is a repulsive octet state, 
one can then wonder if it would not be energetically 
favourable to consider a P-wave for the qq instead of 
exciting the gluon. The trial function respecting this as- 
sumption is (ps{p,^) in (|24|) . This ansatz (l q q = 1), 
not investigated in Ref. JlCJ, corresponds to an electric 
gluon hybrid with quantum numbers J , J G [0,3] or 
J-+.JG [0,2]. 

First, we used only one single trial function depending 
on the hyperradius. We considered the same Hamiltonian 
H with the same value for tr, a s , m c , /i and fi g . We found 
a mass Mq = 4.225 GeV which is already lower than 
the mass of a magnetic gluon hybrid. This time, the 
variational parameter has the value (3 — 0.633 GeV 1 / 2 . 

As a next step, we diagonalized with one 

Gaussian trial function with two parameters. We found a 
mass Mg^ — 4.137 GeV with only one Gaussian function 
ifB^Pi A). The optimal values for the parameters read: 
ax = 0.119 GeV and b x = 0.232 GeV. As in the case of a 
magnetic gluon, we found a mass roughly 100 MeV lower 
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than with one function with the hyperspherical formal- 
ism. We plot the spatial part of the wave functions in 
Figs 0] and [5] (remind the factor A in the solid spherical 
harmonic) . 

Wave function given by the hypersherical formalism 



0.1 .. 
0.08. 



0.06- 
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FIG. 4: Wave function (|17p for an electric gluon hybrid meson. 



TABLE II: Wave function coefficients for electric gluon 



i 


on icr 3 






1 


38.013 


0.119838 


0.232252 


2 


0.948 


0.051007 


0.087008 


3 


—0.977 


039315 


393460 


4 


5.333 


0.319606 


0.714934 


5 


-18.77 


0.091608 


1.320736 


6 


-5.593 


0.398742 


0.209574 


7 


9.742 


0.070955 


1.500372 


8 


55.451 


0.257114 


0.710774 


9 


2.330 


0.066225 


0.219531 


10 


34.224 


1.050299 


2.419412 



the simple Gaussian and the true is five time lower than 
the difference between the two first approximations. The 
wave functions for N — 10 is shown in Fig. [6] The rel- 
ative error between the first trial function (fTT|) and the 
eigenfunction is 29% and the relative error with the single 
Gaussian approximation is 26%. 

Wave function with N = 10 functions in the basis 



Wave function with N = 1 function in the basis 




FIG. 5: Wave function (|23[) for an electric gluon hybrid meson. 

The next step is obviously to check how accurate is 
the one Gaussian approximation by looking for the exact 
eigenvalue of -ff (/i, /i 9 ) (by exact, we mean a value stable 
at 1 MeV). When increasing N, the dimension of our 
basis, we reach a stable mass from N — 10. We check up 
to N = 30 the stability of our eigenvalue. The mass of 
electric gluon hybrid reads 

M { b 0) = 4.115 GeV. (28) 

The coefficients and parameters of each functions can 
be read in Table HT1 Once again, the difference between 



We then conclude that the lowest exotic hybrid mesons 
, 1 % 3 h arc dominated by the electric gluon com- 
ponent. 

C. P-wave hybrid meson 

We have identified two main components in the low- 
est P-wave hybrid mesons. The electric gluon compo- 
nent was resolved to be lighter that the magnetic gluon 
component since the quarks are in an octet state. The 
physical wave function should be a mixing of those two 
situations since our Hamiltonian couples the two config- 
urations. We expect then a lower mass for the ground 
state since coupling two channels repels each other. The 
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more general ansatz for a P-wave hybrid meson is 

N 

P) = a * ex P ( _a * A2 ~ fei ^ 2 ~ 2CiA ' <°) ^g) 

x yiu(di\ + eip), 

and correspond to h , 1 , 1 h ,2 h . Using this 
ansatz, we indeed find a lowest value for the ground state 
of H(jjL,fj, g ): 

M = 4.068 GeV. (30) 

The resulting mass is indeed lower than in the previous 
cases. 

D. Semi-relativistic Hamiltonian 

We showed in the previous sections that the ansatz for 
the wave function is essential to find the correct eigen- 
values of an operator. One simple exponential of the hy- 
perradius may lead to a mass around 100 MeV above the 
true value. But, in this work, we assume that the correct 
Hamiltonian for the qqg system is H since the einbein are 
operators and not numbers. As a consequence, eigenvalue 
of H(/j,,fx g ) are upper bounds of eigenvalues of H. The 
auxiliary fields were introduced to get rid of the square 
roots in Hq. Treating /i and fi g as ordinary numbers in- 
stead of operator causes an error on the energy around 
5% (23|. Let us note also that, the more einbein there 
are, the less accurate is the approximation. We can avoid 
this approximation by diagonalizing H given by ([I]). 

The Fourier transform of the general trial function (|22p 
has also a Gaussian shape. We can then easily compute 
matrix elements for operator in momentum space such as 
the semi-relativistic kinetic operator [20j . With a basis of 
Gaussian functions, we do not need to introduce auxiliary 
fields. This avoids also to determine the optimal values 
of /i and fig. 

We are in position to find the accurate eigenvalues of 
the semi-relativistic Hamiltonian H (JTJ). We computed 
the masses of the lightest hybrid with magnetic (<pa) and 
electric (pa) ghions with N = 1 and N — 10 functions in 
the basis. The masses are below the ones obtained with 
the auxiliary field approximation and read 

M^2 = 4.402 GeV, M$ = 4.017 GeV, (31) 
M^ 0) = 4.361 GeV, M^ 0) = 3.970 GeV. (32) 

The coefficients and parameters of each functions can be 
read in Table IIVI for a magnetic gluon and in Table [V] 
for an electric gluon. The relative difference between the 
real values of H and H(fx,fj, g ) is less than 4%. With the 
wave functions, we computed the expectation values (for 
N = 20) 

Mo = (*l v /p2 + m2|$), Mff0 = (*|^|*>. (33) 



In both cases, magnetic (A) and electric (B) gluons, the 
expectation values of these operators are compatible with 
the parameters (fT9]) of ficj |: 

Mo = 1-652 GeV, nfo = 0-976 GeV, (34) 
/if = 1.736 GeV, ^ = 0.750 GeV. (35) 

As expected, < /if and /t^ > /j,^ since in the mag- 
netic case, the P-wave is reported on the gluon. 



Wave function for the semi-relativistic Hamiltonian 




FIG. 7: Wave function for the semi-relativistic Hamiltonian 
(|22[) with 10 Gaussian functions ((23} for a magnetic gluon 
hybrid meson. 



Wave function for the semi-relativistic Hamiltonian 
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FIG. 8: Wave function for the semi-relativistic Hamiltonian 
l]22p with 10 Gaussian functions (|23[1 for an electric gluon 
hybrid meson. 

The final step is, of course, to diagonalize H without 
any assumption on the wave function, i.e. with (|29|) . The 
ground state mass decreases once again and we obtain 

M = 3.895 GeV. (36) 
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TABLE III: Spectra of JTOj and (0 with different wave func- 
tions. Parameters (|18[l are used. All masses are in GeV. 





Magn. (A) 


Elec. (B) 


Full 




4.573 


4.225 




H(», N ),N=l 


4.462 


4.137 


4.137 


H(fi,fj, s ), N — > 00 


4.445 


4.115 


4.068 


H,N = 1 


4.402 


4.017 


4.017 


H, N -> 00 


4.361 


3.970 


3.895 



We summarize all results in Table Hill namely the spectra 
of H(n,iJL g ) and H for magnetic (A) and electric (B) 
gluons, but also for a more general P-wave ccg hybrid 
mesons (Full). In each case, we displayed the results 
for TV = 1 trial Gaussian functions and the true values 
obtained typically with N = 10 trial functions. In the 
general P-wave cases, the N = 1 give always the electric 
gluon mass since the latter is the main component of the 
wave function. In view of this Table, we notice a gain of 
600 MeV when relaxing all hypothesis ! 



TABLE IV: Wave function coefficients for magnetic gluon 
with the semi-relativistic Hamiltonian. 



i 


a* KT 2 


Oi 


bi 


1 


3.138 


0.086588 


0.205722 


2 


6.262 


0.216618 


0.781620 


3 


0.179 


0.055116 


0.070758 


4 


-0.287 


0.608667 


0.201197 


5 


108.510 


0.875692 


0.739526 


6 


3.673 


0.161373 


0.426376 


7 


-0.742 


0.059422 


1.098722 


8 


-2.084 


0.339798 


0.350189 


9 


39.117 


1.054584 


0.569019 


10 


-145.020 


0.939776 


0.679936 



TABLE V: Wave function coefficients for electric gluon with 
the semi-relativistic Hamiltonian. 



i 


a t 10" 2 




bi 


1 


3.672 


0.125160 


0.214905 


2 


0.130 


0.059228 


0.058032 


3 


89.281 


0.342251 


0.796071 


4 


3.431 


0.580957 


0.262103 


5 


-0.666 


0.150242 


1.634497 


6 


-8.155 


0.491529 


0.358052 


7 


21.727 


1.401183 


3.854743 


8 


-75.356 


0.341822 


0.825713 


9 


-0.076 


0.085946 


0.747611 


10 


-3.160 


2.586533 


0.989884 



In order to compare qualitatively the different wave 
functions, we display in Fig. [5] the wave functions for a 
magnetic gluon at A = and in Fig.rTUlthe wave functions 
for an electric gluon at p = 0. 

For the sake of completeness and in order to compare 
all the wave functions quantitatively, we add the contour 



lines of all the eight wave functions described in the text. 
Those plots are displayed in Figs. [T^ [H EE OS QZ1 [H 
US and HOI 

Wave function for magnetic gluon at 1=0 



0.07 




1 2 3 4 5 6 

IPl 



FIG. 9: Wave function for a magnetic gluon at A = 0. 



Wave function for electric gluon atp=0 



0.16 




1 2 3 4 5 6 



FIG. 10: Wave function for an electric gluon at p — 0. 



The mass of the hybrid with an electric gluon is re- 
solved to be lighter than the magnetic one. Our simple 
Hamiltonian mixes the two states and the mass of the 
lightest states has clearly a bigger electric component. 
However, the decay modes of these two states are differ- 
ent. An electric gluon allows the hybrid to decay into 
a S-wave D*D*. This may induce difficulty to single 
out a hybrid interpretation with respect to conventional 
charmonium. I would be interesting to investigate more 
deeply the decay properties of those states which could 
guide experimentalist in the search for hybrid charmonia. 



9 



PARAMETER INFLUENCE 



The previous sections emphasized on the fact that var- 
ious approximations may lead to overestimation of hy- 
brid masses. However, the single Gaussian approxima- 
tion was resolved to be an acceptable approximation. 
We now turn our attention to this specific approximation 
and test its evolution for different values of the param- 
eters for the semi-relativistic Hamiltonian ([TJ) . We then 
change one parameter and keep the others to their "op- 
timal" values (fT8|) . The results are displayed in Figs fTTI 
and [12] for respectively the evolution of the gluon mass 
and the string tension. Varying the gluon mass is in- 
teresting since effective approaches used a non-vanishing 
value for the gluon mass in the kinetic energy, see for 
instance Ref. [3]. For wide ranges of the parameters, the 
single Gaussian approximation is resolved always to give 
an overestimation of ~ 40 MeV constant in the intervals. 
For a wide range of variation of a, the single Gaussian ap- 
proximation lie around 20-50 MeV above the true mass. 
For hybrid meson systems, the approximative wave func- 
tion (|23p is then robust under parameter evolution. We 
checked that we obtain similar results for electric gluon 
wave functions. 



True mass 

1 Gaussian approximation 




FIG. 11: Hybrid mass evolution according to gluon mass. 



VI. CONCLUSION 

During the main part of this paper, we were concerned 
by Kalashnikova and Nefediev's model [l(| in which the 
starting point was Hamiltonian (fT]). This model shares 
the constituent model's typical features. The kinetic en- 
ergy is the semi-relativistic expression valid for massless 
particles (the constituent gluon). A linear (in fact a V- 
junction) plus Coulomb term is used for the potential. 
The parameters (charm quark mass, string tension and 




0.17 
Sigma (GeV 2 ) 



FIG. 12: Hybrid mass evolution according to string tension. 



strong coupling) were determined on a similar model for 
charmonium. The parameters (|18|) reproduce a quarko- 
nium spectrum in agreement with experimental data and 
are expected to predict an acceptable value for hybrid 
masses. 

However, solving a three-body systems may require hy- 
pothesis: On one side in the Hamiltonian by introducing 
auxiliary fields to get of rid of square roots; on the other 
side on the wave function by imposing particular val- 
ues for internal quantum numbers and/or restricting the 
space spanned by the wave function to a restricted sub- 
space of the total Hilbert space. 

We demonstrated in this work that one should be 
very careful with the approximations made to solve the 
eigenvalue problems numerically. Indeed, auxiliary fields, 
hyperspherical formalism with only one function, Fock 
space reduction (magnetic/electric gluons) may causes 
an overestimation of the mass by amount of 100-600 
MeV. However, in each cases of study (magnetic or elec- 
tric gluon), the less worst approximation was the only 
Gaussian function. For the same Hamiltonian, the sin- 
gle Gaussian overestimated the mass by only ~ 30 MeV 
where the hyperspherical formalism with only one func- 
tion gave always a discrepancy around 130 MeV. 

In their conclusion, the authors of fiol ] pointed out 
a numerical vector hybrid meson mass of 4.397 GeV, 
substantially higher than the experimental candidate 
V(4260). The improvements of the numerical method 
developed in the present paper clearly favor the hy- 
brid meson interpretation of the candidate. Indeed, the 
real eigenvalue of the Hamiltonian PH|) used in [lfj is 
130 MeV below the approximation with the hyperspher- 
ical formalism, the exact discrepancy with the vector 
X(4260). We stressed, nevertheless, that the correction 
to the bare mass were computed with the approxima- 
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tive wave function and should slightly differ. It would 
be interesting to investigate how strong is the difference. 
Moreover, we have shown that the mass remains more or 
less stable under the approximations considered in this 
paper but we think that the decays properties, depend- 
ing strongly of the wave functions, should change more 
that the mass. 

It is also worth mentioning the two other vector states 
Y(4325) from BaBar [|| and Y(4360) from Belle [H[. 
Their mass lie close to the X(4260) and mixing may cause 
a mass shift with respect to pure ccg hybrid states. How- 
ever, the purpose of this paper was the investigation of 



the numerical procedure for a simple Hamiltonian for 
pure ccg. A detail investigation of the possible mixing 
with other states would require a more complicated de- 
scription of the system and is beyond this study. 
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FIG. 13: Wave function (|17[1 for a magnetic gluon hybrid 
meson. 



FIG. 14: Wave function (|23[1 for a magnetic gluon hybrid 
meson. 




FIG. 16: Wave function for the semi-relativistic Hamiltonian 
(1221) with 10 Gaussian functions (|23|l for a magnetic gluon 
hybrid meson. 



FIG. 17: Wave function (|17p for an electric gluon hybrid me- 
son. 




FIG. 18: Wave function (|23|) for an electric gluon hybrid me- 
son. 




FIG. 19: Wave function with 10 Gaussian functions 
for an electric gluon hybrid meson. 
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FIG. 20: Wave function for the semi-relativistic Hamiltonian 
(|22[) with 10 Gaussian functions (|23|l for an electric gluon 
hybrid meson. 



